function  nest(d, c, x)
  implicit none
  real :: nest
  real, dimension(d+1), intent(in) :: c
  real, intent(in) :: x
  integer, intent(in) :: d
  nest = 0.0
end function nest

  

program test_nest
  implicit none
  real :: nest
  real, dimension(2) :: c = (/1.0,2.0/)
  integer :: d = 1
  write(*,*) nest(d, c, 0.5)
end program test_nest
